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We investigate how the infrared behavior of electric and magnetic gluon propagators in Landau 
gauge is affected by temperature. More precisely, we perform large-lattice simulations in pure 
SU(2) gauge theory around the transition temperature Tc and study the longitudinal (electric) 
and transverse (magnetic) gluon propagators in momentum space, proposing the calculation of 
screening masses through an Ansatz from the zero-temperature case. Going from zero to nonzero 
temperature, we see that the longitudinal gluon propagator is enhanced, with an apparent 

plateau value in the infrared, while the transverse propagator Dt{p) gets progressively more 
infrared-suppressed, with a clear turnover in momentum at all nonzero temperatures considered. 
Our data allow us to associate what was previously seen as a peak in the infrared value of Di{p) 
at Tc to severe finite-size effects along the temperature direction. In particular, a temporal lattice 
extent Nt seems to be needed to study the electric sector around the transition. Once these 
systematic errors are eliminated, the infrared behavior of the longitudinal propagator appears to 
be rather independent of the temperature below the transition. Above Tc, the infrared value of 
Dl{p) starts to decrease monotonically with the temperature. 



The XXVIII International Symposium on Lattice Field Theory, Lattice2010 
June 14-19, 2010 
Villasimius, Italy 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Massive Landau-Gauge Propagators at Finite Temperature 



Tereza Mendes 



1. Introduction 

At high temperatures, deconfinement is expected to be felt in the long-distance behavior of 
correlation functions — such as the (real-space) longitudinal gluon propagator — as an exponential 
falloff with the distance, defining a screening length and conversely a screening mass [|l]]. Although 
this predicted behavior has been established at high temperatures down to around twice the critical 
temperature |^, it is not clear how a screening mass would show up around Tc. At the same 
time, recent studies of zero-temperature Landau-gauge propagators on large lattices have shown a 
(dynamical) gluon mass at T = (see e.g. for a discussion), in agreement with the so-called 
massive solution of Schwinger-Dyson equations iQ]. One can try to use this knowledge to define 
temperature-dependent masses for the region around the critical temperature. In the following, we 
review briefly the lattice results for the zero-temperature case, discuss the expected behavior for 
nonzero temperature and show our preliminary results for the gluon propagator on large lattices for 
several values of the temperature, drawing our conclusions. A more detailed analysis and additional 
data will be presented elsewhere [Q]. 

2. Expected Behavior 

At zero temperature. Landau-gauge gluon and ghost propagators are expected to behave ac- 
cording to the so-called Gribov-Zwanziger scenario, which is based on restricting the gauge con- 
figurations to the region delimited by the first Gribov horizon, where the smallest nonzero eigen- 
value of the Faddeev-Popov matrix goes to zero ||^. Essentially, the infinite- volume limit favors 
configurations on the first Gribov horizon and, as a result, the ghost propagator should become 
infrared-enhanced, inducing long-range effects in the theory. (In Coulomb gauge, the restriction to 
the first Gribov region causes the appearance of a confining color-Coulomb potential.) In this sce- 
nario, formulated for momentum-space propagators, the long-range features needed to explain the 
color-confinement mechanism are thus manifest in the ghost propagator, whereas the momentum- 
space gluon propagator D{p) is suppressed in the infrared limit. Such a suppression is associated 
with violation of spectral positivity, which is commonly regarded as an indication of gluon con- 
finement. In fact, Z)(0) is originally expected to be zero, corresponding to maximal violation of 
spectral positivity. The parametrization of this behavior as a propagator with a pair of poles with 
conjugate complex masses was proposed by Gribov, in connection with his study of gauge copies. 

Lattice studies have confirmed the suppression of the gluon propagator in the infrared limit and 
the enhancement of the ghost propagator at intermediate momenta. However, once the investigated 
lattice sizes were large enough, it became clear that the standard procedure for gauge fixing and 
simulations is not compatible with the original scenario in the deep infrared regime. Indeed, the 
gluon propagator attains a finite value as the momentum is taken to zero and the enhancement of the 
ghost propagator does not persist in this limit. We note the very large lattice sizes employed in order 
to observe such a behavior, L 20 fm and larger. The status of these zero-temperature simulations 
has been recently reviewed in |^]. A good fit of the massive behavior for the gluon propagator is 
obtained from the Gribov-Stingl form, which generalizes the Gribov form described above (see e.g. 
[Q]). In any case, violation of reflection positivity for the real-space gluon propagator is observed 
for all cases studied |^. 
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As the temperature T is turned on, one expects to observe Debye screening of the color charge, 
signaled by screening masses/lengths that can in principle be obtained from the gluon propagator 
[10]. More specifically, chromoelectric (resp. chromomagnetic) screening will be related to the lon- 
gitudinal (resp. transverse) gluon propagator computed at momenta with null temporal component, 
i.e. with /?o = (soft modes). At high temperatures, we expect the real-space longitudinal propa- 
gator to fall off exponentially at long distances, defining a (real) electric screening mass, which can 
be calculated perturbatively to leading order. Also, according to the 3d adjoint-Higgs picture for 
dimensional reduction, we expect the transverse propagator to show a confining behavior at finite 
temperature, in association with a nontrivial magnetic mass (see e.g. [^). The ghost propagator, 
on the other hand, should not depend on T . We note that these propagators are gauge-dependent 
quantities, and the (perturbative) prediction that the propagator poles might be gauge-independent 
must be checked, by considering different gauges. 

The above expectations have been checked and confirmed for the gluon propagator at high 



T for various gauges [□, 11]. The behavior of Landau-gauge gluon and ghost propagators around 



the critical temperature Tc has been investigated in [12]. The study showed a stronger infrared 
suppression for the transverse propagator than for the longitudinal one, confirming the dimensional- 
reduction picture also at smaller temperatures. [We note here that a very recent study [ [131 ] discusses 
whether this suppression is consistent with Dt (0) = and investigates Gribov-copy effects for the 
propagators.] It was also found that the ghost propagator is insensitive to the temperature, as 
predicted. For the longitudinal gluon propagator, a very interesting behavior was seen: the data 
approach a plateau (as a function of the momentum) in the infrared region and, as a function of 
temperature, this plateau shows a sharp peak around the critical temperature. The exact behavior 
around T^- (e.g. whether the peak turns into a divergence at infinite volume) could not be determined, 
since relatively small lattices were used. All studies mentioned so far are for SU(2) gauge theory. 
The momentum-space expressions for the transverse and longitudinal gluon propagators Dt{p) 
and Di{p) can be found e.g. in [p^]. 

Recently, in [14], further simulations around Tc confirmed the above results, and lattice data 
for the gluon propagator were used to construct an order parameter for the chiral/deconfinement 
transition. More precisely, the authors use a much finer resolution around Tc and consider the 
SU(2) and SU(3) cases. A check of their calculation is done for the electric screening mass, taken 
as Di{0)^^l^ and extracted from the data, where only the p = raw data point is used. The 
considered lattice sizes are still moderate. 

Of course, even if an exponential fit to the longitudinal gluon propagator works at high tem- 
perature, it is not obvious that this should hold at T > Tc. One should therefore consider more 
general fits. At T = 0, the momentum-space propagator is well fitted by a Gribov-Stingl form (see 
e.g. allowing for complex-conjugate poles 

This expression corresponds to two poles, at masses = a ± ib, where m = niR + imj. The 
mass m thus depends only on a, ^ and not on the normalization C. The parameter rj should be 1 
if the fitting form also describes the large momenta region (from our infrared data we get Tj ^ 1). 
For consistency with the usual definition of electric screening mass, we expect to observe mj — > 
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(Z? — )■ 0) for the longitudinal gluon propagator at high temperature. Clearly, if the propagator has the 
above form, then the screening mass defined by Dl{Q)-^I'^ = ^/{a^ + b'^)/C mixes the complex 
and imaginary masses niR and mj and depends on the (a priori arbitrary) normalization C. 

3. Results 

We have considered the pure SU(2) case, with a standard Wilson action. For our runs we 
employ a cold start, performing a projection on positive Polyakov loop configurations. Also, gauge 
fixing is done using stochastic overrelaxation and the gluon dressing functions are normalized to 1 
at 2 GeV. We take j3 values in the scaling region and lattice sizes ranging from Ns = 48 to 192 and 
from A^, = 2 to 16 lattice points, respectively along the spatial and along the temporal directions. 
The temperature is given by T = 1 /Nt a. 

All our data have been fitted to a Gribov-Stingl behavior, as described in the previous section 



(see Eq. [2. 1[ ). These fits are shown below in all plots, whereas a detailed discussion of the associated 
masses niR, mj will be presented elsewhere [||]. We generally find good fits to the Gribov-Stingl 
form (including the full range of momenta), with nonzero real and imaginary parts of the pole 
masses in all cases. For the transverse propagator Dt{p), the masses m« and ni/ are of comparable 
size. The same holds for Di{p), but in this case the relative size of the imaginary mass seems 
to decrease with increasing temperature. We also looked at the real-space propagators, finding 
clear violation of reflection positivity for the transverse propagator at all temperatures. For the 
longitudinal propagator, positivity violation is observed unequivocally only at zero temperature 
and for a few cases around the critical region, in association with the severe systematic errors 
discussed below. 

Our runs were initially planned under the assumption that a temporal extent Nt = 4 might be 
sufficient to observe the infrared behavior of the propagators. (Our goal was, then, to increase Ns 
significantly, to check for finite-size effects.) For this value of A^f, the chosen jS values: 2.2^12, 
2.299, 2.313, 2.333, 2.505796 yield temperatures respectively of 0.968, 1.0, 1.04, 1.1, 1.936 times 
the critical temperature Tc. As is clear from Fig. [l] below, the assumption that Nt = 4 might be 
enough is not verified for the longitudinal propagator around the critical temperature, especially in 
the case of larger Ns- 

Indeed, as Ng is doubled from 48 to 96 and then to 192, we see that the infrared value of 
Dl{p) changes drastically, resulting in a qualitatively different curve at Ng = 192, apparently with 
a turnover in momentum. Also, in this case the real-space longitudinal propagator manifestly 
violates reflection positivity. We took this as an indication that our choice of A^f = 4 was not valid 
and therefore considered larger values of N, obtaining a clearer picture of the critical behavior of 
Di{p). As seen in Fig. |l] (top), once we use large enough values of Nt, the curve stabilizes (within 
statistical errors) for four different combinations of parameters. In particular, the two curves at fixed 
physical volume (the yellow and the gray curves) agree very nicely. Note also that the (orange) 
curve corresponding to the smallest physical spatial size (i.e. 4 fm), may show mild finite-size 
effects. As seen in the bottom plot, the finite-physical-size effects are more pronounced for the 
transverse propagator, which does not seem to suffer from the same small-A'^r effects. In this case, 
we see clearly the strong infrared suppression of the propagator, with a turnover at around 400 
MeV. 
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Figure 1: Longitudinal (top) and transverse (bottom) gluon propagator at T^, for various lattice sizes and 
values of j3. Values for x Nt, j3, lattice spacing a and spatial lattice size L (both in fm, in parentheses) 
are given in the plot labels. The resulting temperature is about 302 MeV. 
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Figure 2: Longitudinal and transverse gluon propagators at T = (left) and T = Tc/2 (right). Values for 
A^,-' X A';, j3, lattice spacing a and spatial lattice size L (both in fm, in parentheses) are given in the plot labels. 
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Figure 3: Longitudinal gluon propagator for T < Tc (left) and T > Tc (right). Values for A^^^ xNt, j5, lattice 
spacing a and spatial lattice size L (both in fm, in parentheses) are given in the plot labels. 



The above systematic errors were not observed for other simulated values of the temperature 
(except for the data at 0.968 Tc, not shown here), including the temperature values just above Tc. 
We also did runs at T = and T ?a Tc/2, for which we show combined Di{p) and Dt{p) data in 
Fig. g We see a jump in Dl{p) as temperature is switched on, while Dt{p) decreases slightly, 
showing a clear turnover point at around 350 MeV. In Fig. ^ we collect data for Dl{p) at several 
temperatures (for T < Tc, for clarity, we show only valid lattices with the largest physical size). 
Note that the curve remains unchanged (within errors) from Tc/2 to Tc. Above Tc, there is a steady 
drop. For all values of T, Di{p) seems to reach a plateau at small p. 
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4. Conclusions 

The transverse gluon propagator Dt{p) shows infrared suppression and a turnover in momen- 
tum (in agreement with the dimensional-reduction picture) at all nonzero temperatures considered. 
The longitudinal propagator Di{p), on the contrary, appears to reach a plateau at small momenta. 
The data for Di{p) are subject to severe finite-A/; effects at T « r^. As a result, only lattices with 
Nt >8 seem to be free from systematic errors. After these errors are removed, we see an infrared 
value about 50% smaller than before. This suggests that there is no jump in the infrared value of 
Dl{p) as T ^ Tc from below and that the drop after is significantly smaller than previously 
observed. [We note that all previous studies of Dl{p) around have employed Nt < 4.] Therefore, 
the qualitative behavior of Di{p) around the transition has to be revised. 

We have obtained good fits of our data to a Gribov-Stingl form, with comparable real and 
imaginary parts of the pole masses, also in the longitudinal-propagator case. This is in contrast 
with an electric screening mass defined by the expression Di{0)^^l^, which moreover may contain 
significant finite-size effects. In that respect, we plan to consider the upper and lower bounds for 



£)(0) introduced in [15 1, to investigate the infinite- volume limit of the gluon propagator at nonzero 
temperature. 
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